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ABSTRACT 

In a former paper (Garcia- Rissmann ct al. 2005; hereafter Paper I), we have presented 
spectra of 64 active, 9 normal and 5 Starburst galaxies in the region around the near- 
IR Calcium triplet absorption lines (CaT) and the [SIIIJA9069 line. In the present 
paper we analyze the CaT strength (WcaT)i and kinematical products derived in that 
study, namely stellar (a^) and ionized gas (cgas) velocity dispersions. Our main results 
may be summarized as follows: (1) Seyfert 2s show no sign of dilution in VFcaT with 
respect to the values spanned by normal galaxies, even when optical absorption lines 
such as the Call K band at 3933 A are much weaker than in old, bulge-like stellar 
populations. (2) The location of Seyfert 2s in the WcaT-WcaK plane is consistent 
with evolutionary synthesis models. The implication is that the source responsible 
for the dilution of optical lines in these AGN is a young stellar population, rather 
than an AGN featureless continuum, confirming the conclusion of the pioneer study 
of Terlevich, Diaz & Terlevich. (3) In Seyfert Is, both iy[siii] and M^caT tend to 
be diluted due to the presence of a non-stellar component, in agreement with the 
unification paradigm. (4) A comparison of cr+ with (jgas (obtained from the core of 
the [SIII] emitting line) confirms the existence of a correlation between the typical 
velocities of stars and clouds of the Narrow Line Region. The strength and scatter 
around this correlation are similar to those previously obtained from the [OIIIJA5007 
line width. 



1 INTRODUCTION 



The stellar absorption lines from the near-IR Calcium triplet 
(CaT) made their debut in the field of active galactic nu- 
clei (AGN) with the work of Terlevich, Di'az & Terlevich 
(1990, hereinafter TDT). At that time, it was thought that 
the optical spectrum of Seyfert 2s contained a featureless 
continuum (FC) from the AGN, which, in light of unifica- 
tion scenarios (Antonucci 1993), should be associated with 
scattered light from the hidden Seyfert 1 nucleus (Cid Fer- 
nandes & Terlevich 1995). TDT realized that this idea could 
not be fully correct, since their data indicated that, unlike 
for absorption lines in the optical, the strength of the CaT 
in these sources shows no signs of dilution with respect to 
the values found in normal galaxies. Their interpretation of 
this fact was that what was called a "non-stellar FC" was 
in fact a stellar component associated to young stellar pop- 



ulations. This conclusion was thoroughly confirmed by nu- 
merous studies in the past decade (see Cid Fernandes et al. 
2004 for a review). 

Since then, interest in the CaT shifted towards its use as 
a tracer of stellar kinematics, particularly in AGN. The loca- 
tion of the CaT in a relatively clean spectral region makes it 
an ideal feature to measure stellar velocity dispersions ((t*), 
which trace the gravitational potential of the host galaxy's 
bulge. In an influential work. Nelson & Whittle (1996; NW) 
used the CaT in a comparison between stellar and gaseous 
kinematics of AGN, finding that the typical velocity Nar- 
row Line Region (NLR) clouds (cTgas) correlates with cr*, 
albeit with significant scatter. More recently, a strong cor- 
relation between cr* and black-hole mass (Tremaine et al. 
2002 and references therein) indirectly boosted interest in 
observations of the CaT (Nelson 2000; Botte etal. 2004). 

For completely different reasons, recent work on nor- 



2 Vega etal. 



mal galaxies has also raised the interest in the CaT. Empir- 
ical investigations have revealed rather puzzling behavior of 
the CaT in bulges and elliptical galaxies. Firstly, detailed 
population synthesis models tend to over-predict the CaT 
strength (Saglia etal. 2002; Cenarro etal. 2003, 2004), par- 
ticularly for giant ellipticals, though for dwarf ellipticals the 
match between data and models is satisfactory (Michielsen 
et al. 2007) . Secondly, while classic metallicity tracers like 
the Mg2 index are known to correlate with a* (eg, Terlevich 
etal. 1981), tracing the well-known mass- metallicity rela- 
tion, the strength of the CaT appears to be anti- correlated 
with (7* (Cenarro et al. 2003, 2004; Falcon-Barroso et al. 
2003; Michielsen etal. 2003). The interpretation of these re- 
sults is still not clear, and the situation is likely to be even 
more complex in systems with varied star formation histo- 
ries like AGN (Cid Fernandes etal. 2004; Wild etal. 2007). 

One thus sees that, more than 15 years after TDT in- 
troduced the CaT in AGN research, there are still plenty of 
reasons to keep studying it. With this general motivation, in 
Paper I we have carried out a spectroscopic survey of Seyfert 
galaxies in the region including the CaT and the [SIIIJA9069 
emission line. Paper I concentrated in the presentation of the 
data and the derivation of four main data products: stellar 
velocity dispersions (cr*), [SIII] emission line widths ((T[siii], 
which is representative of the highly ionized component of 
the gas of the Narrow Line Region in AGNs), [SIII] equiv- 
alent widths (VK[siii]), and the CaT strength {Wcarr)- Here 
we use these data to address the following issues: (1) ex- 
amine the contribution of a non-stellar component in the 
near-IR spectra of Seyfert galaxies and test the consistency 
of VK[siii] and VKcaT data with the unified model, (2) study 
the sensitivity of IVcaT to stellar population properties, (3) 
investigate the connection between NLR and stellar kine- 
matics and CaT strength, and (4) report results of spatially 
resolved spectroscopy, not reported in Paper I. 

To study these questions, we complement the data in 
Paper I with values of the equivalent width of the Call K 
line (W^CaK) obtained from previous studies. As shown by 
Cid Fernandes etal. (2001), the Call K line is a powerful 
tracer of stellar populations even in type 2 Seyferts. By com- 
paring WcaTT with WcnK we can assess whether the CaT is a 
useful stellar population tracer, and determine whether the 
combination of CaK and CaT strengths is explained as due 
to stars alone or if an FC component is necessary. We also 
add in literature information on the width ((T[oiii]) of the 
[OIIIJA5007 line, which is useful to test whether our results 
are somehow affected by the choice of [SIII]A9069 as a tracer 
of NLR motions. 

This paper is organized as follows. Section [2] presents a 
series of studies related to the CaT strength. After showing 
CaT observational properties in l2.1l we compare WcaT with 
WcaK in 12.21 and investigate whether the location of the 
points in the PFcaT versus WcaK diagram is consistent with 
the existence of an FC at near-IR wavelengths. In|3]we use 
evolutionary synthesis models to track the behavior of both 
WcaT and WcaK as a function of age and metallicity, and 
overlay the models onto the data in the WcaT versus WcaK 
diagram. In|4]we study the spatial behavior of Gi, and WcaT 
for about half of our sample. In Section [S] we analyze stellar 
and ionized gas kinematics: stellar velocity dispersions and 
its relation to CaT are studied in l5.ll while in l5.2l we discuss 



the link between stellar and gas kinematics. Finally, section 
[6] summarizes our main results. 



2 THE CAT STRENGTH: EFFECTS OF 
NON-STELLAR LIGHT 

One of our goals is to evaluate the usefulness of the CaT as a 
stellar population diagnostic. Since our sample is dominated 
by AGN, before discussing stellar populations we must first 
examine to which extent non-stellar light affects our CaT 
measurements. 

As explained in Paper I, we adopt the W^caT and H^caT* 
definitions of Cenarro et al (2001a). They offer two defini- 
tions of the CaT equivalent width: 'CaT' (which we call 
Wcarr), which consists of a sum of the equivalent widths of 
all three CaT lines, and 'CaT*' (called WcaT* here), which 
corrects VFcaT for contamination by Paschen line absorp- 
tion. These equivalent widths are measured with respect to 
a continuum defined by fitting the spectrum in 5 windows in 
the 8474-8792 A range. Because of the presence of unwanted 
features in the observed spectra, we opted to measure WcbTT 
and Wca.T* on the synthetic spectra. On the other hand, our 
spectral base used to fit the observed spectra does not have 
stars with the spectral types with Paschen lines, so we use 
WcaT as equivalent widths measurements for our spectra. 

All nuclear CaT data refers to appertures of ~ 2" x 2" 
(Paper I). To take into account the different velocity disper- 
sions of our galaxies, we broadened each spectrum to match 
the largest value of ct*^ -f cTinst^ of our sample, correspond- 
ing to NGC 3115 {a-, = 275 and ain.t = 56 km/s). This 
correction changes the value of VKcaT by typically 10% with 
respect to the uncorrected values given in Paper I. VKcaK 
measurements from Cid Fernandes et al (2004) are also used 
in our analysis. 

2.1 Seyfert Is x Seyfert 2s 

It is well known that a non stellar light is present in Seyfert 
Is. The contribution of such a component in Seyfert 2s has 
been a topic of much controversy in the past. One way to 
address this issue is to compare Seyfert Is with Seyfert 2s 
through CaT. 

Fig. [T] (top panel) shows the distributions of Wcarr 
for Seyfert galaxies in our sample, with Seyfert Is (includ- 
ing subtypes 1-1.9) marked by the filled areas. Almost all 
Seyfert 2 lie between 4 and 8 A, with a median of 6.5 A. 
M^CaT tends to be smaller in type 1 nuclei, with a median of 
4.8 A. Emission lines from the NLR follow the same pattern, 
with smaller equivalent widths in Seyfert Is than in Seyfert 
2s, as shown for [SIII]A9069 in the bottom panel of Fig. [T] 

Another way of illustrating this result is shown in Fig. 
(2] where we plot the CaT strength as a function the frac- 
tional contribution of power-laws to the continuum in the 
CaT region. This fraction, which we denote xpl, is deduced 
from the spectral fits presented in Paper I, where, besides 
observed stars, power-laws of several slopes were included to 
emulate the effects of a template mismatch in the CJ^, mea- 
surements. Table [1] list xpL, along with the new (corrected) 
VKcaT values mentioned at the beginning of this section and 
their errors. 

The plot shows that xpL is smaller for Seyfert 2s than 
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Name 


Type 


Wc^T (A) 




(%) 


Name 


Type 


W^CaT (A) 


XPL (%) 


NGC 


205 


Normal 


5.68±0.46 




32 


NGC 7410 


Sy 2 


7.10±0.46 


22 


NGC 


526A 


Sy 1.9 


4.46±0.45 




49 


NGC 7469 


Sy 1.5 


2.72±0.45 


65 


NGC 


526B 


Normal 


6.80±0.47 




19 


NGC 7496 


Sy 2 


5.45±0.47 


37 


NGC 


1019 


Sy 1.5 


5.92±0.44 




35 


NGC 7582 


Sy 2 


5.22±0.47 


14 


NGC 


1068 


Sy 2 


5.84±0.48 




24 


NGC 7590 


Sy 2 


7.12±0.43 


17 


NGC 


1125 


Sy 2 


6.89±0.46 




21 


NGC 7714 


Sy 3 


4. 25 ±0.46 


53 


NGC 


1140 


Starburst 


5.46±0.44 




37 


MARK 0001 


Sy 2 


5.64±0.46 


25 


NGC 


1142 


Sy 2 


7.99±0.46 







MARK 0003 


Sy 2 


3.97±0.48 


47 


NGC 


1241 


Sy 2 


7.86±0.43 




18 


MARK 0040 


Sy 1 


4.24±0.46 


38 


NGC 


1365 


Sy 1.8 


0.93±0.46 




85 


MARK 0078 


Sy 2 


6.90±0.45 


15 


NGC 


1380 


Normal 


7.64±0.46 




5 


MARK 0079 


Sy 1.2 


2.64±0.46 


55 


NGC 


1386 


Sy 2 


^ A^ \ f\ A 

7.47±0.46 




5 


MARK 0273 


Sy 2 


7.30±0.44 


10 


NGC 


1433 


Normal 


6.79±0.44 




4 


MARK 0348 


Sy 2 


5.87±0.44 


21 


NGC 


1672 


Starburst 


7.1d±0.44 




14 


MARK 0372 


Sy 1.5 


D.08±0.46 


6 


NGC 


1808 


Starburst 


d.83±0.46 




21 


MARK 0461 


Sy 2 


r' ^ ^ \ r\ AT' 

5.57±0.46 


25 


NGC 


2110 


Sy 2 


6.34±0.47 




19 


MARK 0516 


Sy 1.8 


6.63±0.46 


28 


NGC 


2639 


Sy 1.9 


C ^"7 A \ r\ A A 

o.74±0.44 




27 


MARK 0573 


Sy 2 


7.2d±0.44 


9 


NGC 


2997 


Normal 


7.39±0.45 




16 


MARK 0705 


Sy 1.2 


4.61±0.46 


41 


NGC 


3081 


Sy 2 


6.89±0.46 




24 


MARK 0915 


Sy 1.8 


A f^Ci \ {\ A^ 

4.78±0.47 


28 


NGC 


3115 


Normal 


^7 {\A _i_n A C 

7.04±U.4o 




15 


MAKK lUoD 


by 2 


5.31±0.45 


31 


NGC 


3256 


Starburst 


3.46±0.46 




55 


MARK 1073 


Sy 2 


5.50±0.46 


25 


NGC 


3281 


Sy 2 


6.85±0.45 




19 


ATA T> Ty 1 1 

MARK 1210 


Sy 2 


\ f\ Act 

6.22±0.48 


38 


NGC 


3783 


Sy 1.5 


2.72±0.45 




66 


MARK 1239 


Sy 1.5 


0.92±0.47 


81 


NGC 


4339 


Normal 


c c\-\ \ r\ AC 

6.91±0.46 




6 


KbU 3d2-G08 


Sy 2 


6.87±0.45 


22 


NGC 


4507 


Sy 2 


6.56±0.45 




24 


ESQ 362-G18 


Sy 1.5 


6.14±0.45 


28 


NGC 


4593 


Sy 1 


3.29±0.45 




60 


IC 2560 


Sy 2 


7.32±0.45 


14 




A '7 A Q 


Sy 1 


3.06±0.44 




63 


IC 3639 


Sy 2 


5.80±0.46 


31 


NGC 


4968 


Sy 2 


6.46±0.45 




20 


IC 5169 


Sy 2 


6.88±0.45 


19 


NGC 


5135 


Sy 2 


5.60±0.46 




37 


IRAS 01475-0740 


Sy 2 


5.48±0.47 


36 


NGC 


5929 


Sv 2 


6.08±0.46 




17 


IRAS 04502-0317 


Sy 2 


6.37±0.47 


19 


NGC 


6300 


Sy 2 


7.69±0.44 




17 


MCG -01-24^012 


Sy 2 


6.28±0.47 


27 


NGC 


6814 


Sy 1.5 


3.52±0.47 




47 


MCG -02-08-039 


Sy 2 


7.37±0.46 


5 


NGC 


6860 


Sy 1.5 


5.16±0.43 




48 


MCG -06-.30-015 


Sy 1.5 


4.60±0.45 


47 


NGC 


6907 


Normal 


8.54±0.44 




4 


MCG -1-8-11-11 


Sy 1.5 


0.92±0.44 


84 


NGC 


6951 


Sy 2 


8.28±0.45 




10 


UGC 1395 


Sy 1.9 


5.98±0.45 


33 


NGC 


7130 


Sy 1.9 


6.43±0.45 




25 


UGC 12138 


Sy 1.8 


6.26±0.44 


14 


NGC 


7172 


Sy 2 


6.36±0.44 




25 


UGC 12348 


Sy 2 


6.78±0.45 


12 


NGC 


7184 


Normal 


7.34±0.44 




19 


UGC 3478 


Sy 1.2 


1.28±0.44 


66 


NGC 


7212 


Sy 2 


4.92±0.46 




33 


AKN 564 


Sy 1.8 


3.05±0.49 


32 


ilts of 




and xpL- ' 


Columns lists 


galaxy 


name, 


type of activity, WcaT and 


xpL contribution. All 1 



measured in the same system, as explained in the text. 



for Seyfert Is. Furthermore, the values of both xpl and 2.2 Dilution by an underlying FC: The CaT 
WcaT obtained for Seyfert 2s span approximately the same versus CaK diagram 

range as those found for the non-active galaxies in the sam- 

plg_ While the results reported above suggest one docs not need 

to worry about non-stellar light in the CaT range in Seyfert 
2s, it is desirable to test this idea more conclusively before 
proceeding to an analysis of WcaT in terms of stellar popula- 
tions alone. This can be done combining our CaT data with 
information from a different wavelength range. The equiv- 
alent width of the Call K absorption at 3933 A (Wchk) 
The fact that VFcaT(Seyl) < VKcaT(Sey2) and similarly provides the necessary information. 
W[siii]{Seyl) < W^[s///] (Sey2) is consistent with the unified The CaK line, which originates in old, late-type stars 

model, which predicts a dilution of the equivalent widths from the host's bulge, is very sensitive to dilution by an 
of the absorption and emission lines in sources where the underlying blue continuum. As usual, this continuum can 
nuclear non-stellar continuum is seen directly. The fact that be either due to an AGN or to a young stellar population, 
the range of CaT strengths spanned by Seyfert 2s is similar and the value of WciiK by itself cannot distinguish between 
to that spanned by non-active galaxies indicates that their these two possibilities (e.g., Storchi-Bergmann etal. 2000). 
scattered FC, if present, is too weak to produce a noticeable The combination of WcaK and Wcarr data helps lifting this 
eflfect on the near-IR range. ambiguity. 
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Figure 1. Distributions of the equivalent width of the CaT ab- 
sorption lines (top) and of the [SIIIJA9069 emission line (bottom) 
for Seyfert 1 (filled areas) and Seyfert 2 (empty areas). 
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Figure 2. WcaT versus percentage of power-law contribution to 
the synthesis. Hereafter, symbols are like in this figure. 



The idea, first developed by TDT, is simple: if the CaK 
is diluted by a truly featureless continuum, and if this FC ex- 
tends from the optical to the near IR, then the CaT would be 
diluted too, and by a predictable amount. From observations 
of type 1 Seyferts and quasars, we know that the AGN opti- 
cal continuum is well described by a power-law F^'^ oc v'" , 



with Q in the 1-2 range or, equivalently, oc A~ , with 

/3 = 2 — a. This FC is superposed to the bulge stellar light. 

Suppose that the bulge has intrinsic CaK and CaT 
equivalent widths Wc^k ^^nd WcaTi s-nd that the ratio of 
the stellar continuum fluxes over these two lines is c* = 
'-^^caK/'-^^caT' Trlvlal algebra then leads to the following 
relation between the observed CaK and CaT equivalent 
widths: 



WcaT = WcaT- 



CaK 



(1) 
(2) 



where x^j^ is the fraction of the total continuum at 3933 
A which comes from the FC. 

We modeled a mixing of a FC and bulge continuum in a 
form of "FC -|- SSP" mixing lines, where SSP refers to a Sim- 
ple Stellar Population representative of 10 Gyr old, typical 
of elliptical galaxies and bulges, as obtained from Bruzual 
& Chariot (2003, BC03) models for solar metallicity. Practi- 
cal use of equation [2] requires stipulating fiducial values for 
M^CaK, VK^aT and c\ We adopt (V7^aK, W^Cax) = (21.1, 7.9) 
A and c* = 0.38. 

Equation [2] is over plotted to the data points in the 
M^CaT versus Wca.K diagram (Fig. [3| . The lines indicate the 
bulge + FC mixing lines for different values of the slope 
a and for x^aK decreasing from 100% at (WcaT,W^CaK) = 
(0,0) to 0% at (M^caT,W^CaK) = (V7^aT,W^aaK)- The Wc.K 
measurements come from the optical studies of Cid Fernan- 
des etal. (1998, 2001, 2004). Of the 78 galaxies in our sam- 
ple, 42 overlap with these studies, most of which (31) corre- 
spond to Seyfert 2s. 

We see that some Seyfert nuclei line up along the mix- 
ing lines. However, in order to explain all of the objects, it is 
necessary a contribution of x^^jf of about 50% or greater. 
There are some Seyfert with high VKcaT and low VKcaK that 
are not explained in terms of these mixing lines for a be- 
tween 1 (/3=1) and 2 {/3—0), which are representative of type 

1 AGN. For these objects, only if a were as low as ~ 0.5 
or less the dilution lines would go approximately through 
the observed values, but not even quasars have such an ex- 
tremely blue continuum. Such a blue FC could in principle 
be produced by scattering of the AGN light by dust parti- 
cles. Nevertheless, the required values of the x^^^ fraction 
(of the order of ~ 50%) would imply that scattered broad 
lines in the optical should become visible, in which case the 
galaxy would not be classified as a Seyfert 2 in the first 
place (an argument first put forward by Cid Fernandes & 
Terlevich 1995). Besides that, about half of the Seyfert 2s 
with XQ^jf greater than 50% show High Order Balmer Lines 
(HOBL) in absorption in their spectra (Gonzalez Delgado 
etal. 2001). So, not diluted HOBL argue against a dilution 
by a power-law in these objects. 

We have to point out that "FC -|- SSP" in Fig. [3] are 
computed with a solar SSP values. For another metallicities, 
fidutial numbers are a litUe different (mainly lower H^CaT), 
and the lines would go down for about 1 A. In that way, 
mixing lines in Fig.|3]should be considered as an upper limit 
for "FC + SSP" models. 

We also computed "FC -I- E" models (not shown in Fig. 
[S]), where E refers to NGC 2110, the earliest (t=-3) Seyfert 

2 galaxy in our sample for which we have determined Wca.K 
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5 10 15 

[A] 



25 



Figure 3. Equivalent width of the CaT (from Paper I) plotted 
against the equivalent width of the Call K line (from the litera- 
ture). The lines show the expected dilution curves (equation[2]l for 
bulge plus FC power-laws, with slopes a = 1 and 2 (typical range 
of type 1 AGN), as labeled. Upper scale shows the percentage of 
^CaK contribution. Symbols as in Fig. [2] 



of stellar parameters (Tcff , log g, [Fe/H]) has been addressed 
by several studies (e.g., Diaz, Terlevich & Terlevich 1989), 
and this knowledge has been incorporated into population 
synthesis models which map the CaT dependence on the 
age and metallicity of coeval stellar populations (e.g., Idiart, 
Thevenin & de Freitas Pacheco 1997). In a recent and com- 
prehensive study, Cenarro and coworkers revisited this topic, 
analyzing the CaT in a large library of stars and producing 
state-of-the-art predictions for the behavior of the CaT in 
integrated stellar populations (Cenarro etal. 2001a, b, 2002; 
Vazdekis etal. 2003). 

Despite early hopes that the CaT could serve as a metal- 
licity indicator for galaxies (as it does for globular clus- 
ters, at least below ~ 0.3Zq; Armandroff & Zinn 1988; 
Idiart et al. 1997; Vazdekis et al. 2003) , or as tracers of re- 
cent star-formation through its sensitivity to the presence of 
red super-giants (hereafter RSG; see TDT; Garcia- Vargas, 
MoUa & Bressan 1998; Mayya 1997), observations have re- 
vealed that the CaT shows remarkably little variation in 
strength for galaxies spanning wide ranges in morphological 
and stellar population properties. This small sensitivity is 
qualitatively consistent with the new generation of synthesis 
models (Vazdekis etal. 2003), which show that, for metal- 
licities higher than ~ —0.5 dex, this feature varies little as 
a function of age and metallicity. 

Overall, these results raise serious doubts as to the use- 
fulness of the CaT as a stellar population diagnostic. In this 
section we take another look at this issue, tackling it from 
an empirical perspective. As in the previous section, infor- 
mation on the CaK line will be incorporated in the analysis. 



and WcaT. For this model, (W^^aK. W^Cax) = (13.5,6.34) A, 
and the dilution lines would stay even more distant from the 
data for the same range of a considered before. Not even for 
a as low as ~ (not realistic) the models would along up 
the data. So, our main conclusions about FC -I- bulge mixing 
models for Seyfert 2 nuclei remains valid. 

TDT have performed the same analysis, but using the 
Mgib line at 5175 A instead of CaK (see also Jimenez-Benito 
etal. 2000). They conclude that if the observed dilution of 
MgIb is due to a power-law FC of reasonable slope, then 
the CaT would also be diluted, contrary to their empiri- 
cal result. Here, we considered "FC -I- SSP" mixing models 
with some reasonable combination of Wq^j^, M^caT ^^'^ c* 
values representative of old stellar populations. These mod- 
els shows that dilution by a power-law FC seems to be not 
completely adequate to account for the observed strengths 
of optical and near-infrared ranges simultaneously. Hence, a 
simple mixture of a bulge-like component plus an AGN-like 
FC, envisaged in the early days of AGN research (e.g., Koski 
1978), is just not a viable description of the data. We will 
analyse another explanation for WchK and WcaT behavior 
in the next section. 



3 CAT STRENGTH AS A STELLAR 
POPULATION DIAGNOSTIC 

Now that we considered dilution by non-stellar light in the 
CaT range in type 2 Seyferts, we will proceed to a stellar 
population analysis. 

The behavior of the CaT equivalent width as a function 



3.1 CaT versus CaK 

Let us assume that strongly diluted CaK lines in Seyfert 2s 
are the result of the presence of a Starburst component, in 
line with previous investigations. In particular, Cid Fernan- 
des et al. (2001) propose that Seyfert 2s with WcaK < 10 A 
can be safely identified as composite Starburst -|- AGN sys- 
tems, where the Starburst is unambiguously identified by 
independent features such as the WR bump or PCygni line 
profiles in the UV. Investigating the behavior of WcaT as 
a function of WcaK thus provides a fully empirical test of 
whether the CaT is or is not sensitive to stellar populations. 

A CaT versus CaK diagram was already shown in Fig. 
O where most points correspond to Seyfert 2s. The plot 
shows that WcaT bears little, if any, relation to Wca.K- It 
is therefore hard to evade the conclusion that for the mean 
ages and metallicities of metal rich stellar populations, WchT 
is a poor diagnostic: objects as diverse as Mrk 1210, which 
contains a powerful young Starburst, including Wolf-Rayet 
stars (Storchi-Bergmann, Cid Fernandes & Schmitt 1998), 
and NGC 7172, which is dominated by ~ 10 Gyr stars (Cid 
Fernandes et al 2004) have Wca.T values indistinguishable 
within the errors (6.7 ± 0.4 and 6.9 ± 1.1 A, respectively). 

In Fig. |4]we combine our CaT data with the results of 
the detailed spectral synthesis analysis of type 2 Seyferts 
carried out by Cid Fernandes et al. (2004) for 14 sources 
in common with our sample. We compare their young stars 
(age < 25 Myr) plus "FC" light fractions at 4020 A (xy/fc) 
with both WcaT (left panel), and VKcaK (right). As discussed 
in that paper, which dealt exclusively with spectra in the 
3500-5200 A range, what the synthesis returns as a FC is 
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Figure 4. Percentage of a Starburst component (see text) con- 
tribution to the light in synthetic spectra versus WcaT (left) and 
M^CaK (I'ight) in Seyfert 2s. Note the contrast of the behavior of 
W^CaKi which is very sensititivc to young stars, with the approx- 
imate constancy of WqhT- 
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Figure 5. Evolution of the WcaT*i M^CaT and WcaK absorption 
line indices for SSP, from BC03 models. Long dashed, dashed, 
solid and dotted lines correspond to metallicities of 0.2, 0.4, 1 
and 2.5 Zq, respectively. 



in most cases a reddened Starburst component instead of a 
true non-stellar component, such that Xy/FC is in practice 
a measure of the Starburst strength. 

The near constancy of Wca,T, already evident in Fig. O 
is even more striking in Fig. [l] The inexistent correlation 
between CaT strength and xy/ fc contrasts with the strong 
decrease in WcaK as the Starburst component increases in 
strength. Again, one is lead to the conclusion that, by itself, 
Wca.T is not a stellar population tracer for our galaxies. 

3.2 Evolutionary Synthesis Models in the CaT 
versus CaK plane 

The most natural interpretation of Figs. |3]and|3]is that the 
dilution of optical lines in some Seyfert 2 nuclei is due to a 
young stellar population which causes little or no dilution 
of the CaT. In this section we use models to check whether 
this is a viable scenario. 

To investigate the behavior of the CaT and CaK lines 
as a function of stellar population properties, we have com- 
puted the evolution of VKcaT and WcaK for SSP for 4 metal- 
licities. BC03 evolutionary synthesis models were used for 
this purpose. SSP models are computed with the STELIB 
library (Le Borgne etal. 2003), Chabrier (2003) mass func- 
tion and Padova 1994 tracks (see BC03 for details). Because 
WcaT, WcaT* and WcaK are not provided in the standard 
BC03 distribution, they were computed directly from the 
theoretical spectral energy distributions, as explained at the 
beginning of Section [S] For compatibility with our WchT 
measurements, we broadened all model spectra as explained 
at the beginning of section (2] adopting for the instrumental 



resolution that of the STELIB library (~ 45 km/s at ~ 8600 

A). 

Fig. [5] shows the results for four different metallicities: 
0.2, 0.4, 1 and 2.5 Zq. WcaK is practically null during the 
first ~ 10*'^ yr. After these initial phases, the CaK strength 
increases monotonically with time, which makes it good to 
discriminate very young SSP (younger than ~ 3 x 10* yr) 
from SSP older than ~ 1 Gyr. On the other hand, Wca.T 
starts from small values up to a few Myr, when the first RSG 
appear. The exact timing and duration of this phase depends 
on the metallicity. From ~ 10 Myr to about 1 Gyr, Vl^caT 
oscillates as RSG are replaced by red-giants. Throughout 
these different phases, however, WcaT spans a small range, 
from ~ 6 to 8 A. Taking into account the uncertainties in 
both models and data, this is in practice a constant value. 
Only in the first Myrs of evolution does WcaT have a useful 
diagnostic power, as it assumes much smaller values. Notice 
that in these early phases Wca.T is contaminated by Paschen 
lines, which explains why it does not start from at f = 0. 

In Fig. |6] these models are overlayed to the data in the 
W^CaT versus WcaK diagram. Seyfert Is are not plotted be- 
cause of their non-stellar continuum, which affects both CaK 
and CaT. Numbers along the model tracks mark ages of lO'', 
10'', 10*, 10*^ and 10^° yr. The first thing to notice is that, 
unlike the bulge -I- FC mixing lines in Fig.|3] these models do 
span the region occupied by objects in our sample, showing 
that stellar populations alone are capable of explaining the 
behavior of Seyfert 2s in the CaK versus CaT diagram. Sec- 
ondly, Fig. [6] seems to indicate that sub-solar abundances 
may be required to match the data on Seyfert 2s, partic- 
ularly those with weak CaT. This is somewhat surprising, 
given that we are dealing with the nuclei of early type spi- 
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rals, where Z is expected to be relatively high. One must 
nevertheless recall that these are SSP models, whereas ac- 
tual galaxies contain a mixture of stellar populations. This is 
particularly true for type 2 Seyferts, where previous studies 
have revealed a wide variety of star formation histories (Cid 
Fernandes etal. 2001; 2004; Wild etal. 2007). More realis- 
tic star formation histories can be modeled as a mixture of 
SSP, whose effect on this diagram can be conceived drawing 
imaginary mixing lines connecting any series of points along 
the instantaneous burst tracks. 

3.3 Mixed models: Young plus Old populations 

To have a rough idea of the effects of such a mixture, in Fig. 
Owe repeat the spectral mixing experiment of Fig.[3]but this 
time using a 1 Myr population (to represent ongoing star- 
formation) instead of a non stellar power-law continuum. 
Solar and 2.5 solar BC03 models are used. We further allow 
for the possibility that the young population suffers more 
extinction than the old one, so we consider three cases. Ay 
= 0, 1.5 and 3 mag. Again, we do not plot Seyfert 1 nuclei, as 
they present dilution due to the non-stellar continuum and 
thus are out of the scope of this experiment. This exercise 
is essentially the same as that performed with a power-law 
FC in Fig. [31 except that the continuum shape is not the 
same, and, due to contamination by Paschen lines, WcaT is 
not = for a pure Starburst (see Fig. O middle panel). 

As is evident in the Figure, these mixing lines span the 
observed points for most of the Sy 2 sample. It is important 
to emphasize that only for very young populations one ex- 
pects a relatively low CaT. As shown in the plot, galaxies 
with VKcaT < 6 A and WchK < 10 A require such a popula- 
tion to be adequately modeled. Not even the 0.4 Zq models 
are capable of modeling these objects without a very young 
population. 

It is also important to note that although our 1 Myr + 
10 Gyr mixing models span most of the observed points in 
the CaK versus CaT plane, galaxies with WcaT above ~ 6 A 
can be equally well modeled with an older young population. 
For metallicites > 0.4 Zq, WcaT rises above 6 A as early on 
as ~ 5 Myr, with the appearence of the first RSGs. This also 
happens for continuous star formation models, as shown by 
the solid line in the CaT-CaK diagram. 

Mixing models for > 5 Myr -I- 10 Gyr populations are 
not shown for clarity, but, as can be deduced from Fig. [E] 
they would span most of the upper part of CaK versus CaT 
plane, where many Sy 2s are located. We note in passing 
that TDT originally suggested that the CaT in Seyfert 2s 
signal the presence of RSG in a Starburst. While this may 
well be true in some cases, our modeling shows that there is 
no way of telling it just from the strength of the CaT, since 
we are able to fit the WcaT and WcaK data both with and 
without RSGs. 

Given the simplicity of these models, the overlap be- 
tween data and models in the CaT-CaK diagram is highly 
encouraging. Nevertheless, at the middle of the diagram 
there are two Sy 2 that are far away from our models. The 
two "outliers" are Mrk 3 and NGC 7212. It may be that 
a bona-fide AGN continuum is present in these galaxies, 
which would help explaining their low CaK and specially 
CaT lines. The spectral synthesis analysis of NGC 7212 by 
Cid Fernandes et al. (2004) indicates that ~ 35% of the light 
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Figure 6. Same as Fig. |3] but without Seyfert 1 nuclei and 
over plotting simple stellar populations, as computed from the 
BC03 models. Lines as in Fig. [5] Numbers in stellar lines denotes 
log(Age[yr]). Symbols as in Fig. [2] 



at 4020 A comes from an FC component. They also identify 
a weak broad H/3, strengthening the case for the presence of 
a non-stellar continuum. Furthermore, for both NGC 7212 
and Mrk 3, imaging (Pogge & De Robertis 1993; Kotilainen 
& Ward 1997), spectropolarimetry (Tran 1995), and spec- 
tral modelling (Gonzalez Delgado et al. 2001) all point to 
the existence of an FC component (see Cid Fernandes et al. 
2001 for a detailed discussion). None of these studies find 
compelling evidence for significant on-going star-formation 
in these galaxies. It is therefore likely that their low H^caT 
values (the smallest amongst Seyfert 2s) are indeed due to 
dilution by non-stellar light, explaining why they deviate 
from the region spanned by our purely stellar models in Fig. 

m 

A consistent scenario thus emerges, where the blue con- 
tinuum which dilutes the optical spectrum of some Seyfert 
2s is due to young stellar populations, while their undiluted 
CaT strengths can be explained as due to RSGs or older 
populations. As explained above, in the few cases where the 
CaT is low, a combination of an old plus a very young (pre- 
RSG) population is needed. In these cases, one may expect 
to see Paschen lines in absorption, which is not observed in 
our galaxies, as reported in Paper I. 

To investigate this potential problem we have examined 
the strength of Paschen lines in the theoretical spectra of the 
mixing models shown in Fig. [7] We pay particular attention 
to the Pal4 line at 8598 A, which is strategically placed be- 
tween Ca2 and Ca3. We find that, for Ay = 0, Pal4 only 
becomes visible for > 95% or 85% for Z = 1 and 2.5 

Zq, respectively. These fractions falls to about 75 % is we 
consider 3 magnitudes of differential extinction. These very 
large fractions result from the fact that the young popula- 
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Figure 7. Same as Fig. \E\ now overploting old plus redenned 
young stellar population models for solar and 2.5 solar metallicity. 
Upper scale shows the percentage of the Starburst contribution 
in the CaK region, x^^^^. Mixture lines corresponds to different 
extinction magnitudes and the asterisks denote null (0%, upper 
right) and full (100%, lower left) young contribution. The CSF 
line refers to continuous star formation from ages between 1 Myr 
and 10 Gyr. Dots on this line indicates bursts of 5 and 10 Myr 
old. 



tion is much bluer than the old one, so any mixture where 
young stars contributes significantly to the CaT range, auto- 
maticaly implies that this same component completely dom- 
inates in the optical. Spectral synthesis studies in the optical 
range never find such huge starburst contribution in Seyfert 
2s. Thus, the non detection of Paschen lines may be simply 
due to the fact that old stars dominate the fiux in the near 
IR, even when younger populations dominate the optical 
light. 

A further factor which may be related to the non detec- 
tion of Paschen absorption is filling by an emission compo- 
nent, powered either by the AGN or a circumnuclear star- 
burst. A few of our galaxies do show a hint of Pal4 emis- 
sion, like the Starburst galaxy NGC 7714, where Pal4 has 
about 0.5 A emission equivalent width (see also Gonzalez 
Delgado etal. 1995). In other cases, emission and absorp- 
tion may approximately compensate each other. This seems 
to be the case in the super star clusters of NGC 1569, where 
WR features, high-order Balmer absorption lines and the 
CaT are all clearly detected (Gonzalez Delgado etal. 1997), 
but Paschen lines are not seen either in absorption or emis- 
sion. Hence, even though further studies concerning Paschen 
lines are desirable, their no detection do not invalidate our 
conclusions. 
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Figure 8. Examples of WcaT spatial profiles. Dotted lines mark 
—200 to +200 pc region; dashed lines, the —1 to kpc region. 
Panels (a) and (b) show the dilution of WcaT in Seyfert 1; panels 
(c)-(f) show how flat VFcaT is for Seyfert 2 and Seyfert 2 + 
Starburst. 



4 CAT STRENGTH AS A FUNCTION OF 
DISTANCE TO THE NUCLEUS 

All CaT data reported so far pertains to nuclear regions. For 
34 of the 78 galaxies studied in Paper I our long-slit data 
are of sufficient quality to warrant the extraction of spa- 
tially resolved spectra. Detailed results for all such sources 
are presented in Asari (2006). In this section we summarize 
the main results concerning the spatial behavior of the CaT 
strength. 

Fig.|S]shows W^caT spatial profiles for 6 galaxies, chosen 
to represent Seyfert Is (panels a and b), Seyfert 2s (c and 
d) and composite Starburst -I- Seyfert 2 (e and f ) . There is a 
clear drop in WcaT as one approaches the nucleus in Seyfert 
Is, whereas Seyfert 2s present remarkably fiat Wcajrir) pro- 
files. At a few hundred pc from the nucleus, however, the 
CaT strength becomes similar for both Seyfert types. This 
illustrates once again how sensitive WcaT is to an underly- 
ing FC, in agreement with the results reported in Section 

El 

In a way, these results are arguably more convincing 
proof of the FC effects. Since the analysis in Section [2] was 
based on a comparison of CaT strengths for Seyfert Is and 
2s, or with fiducial values, one could worry that for some 
reason the stellar populations of the two kinds of AGN are 
intrinsically different, invalidating a comparative analysis. 
Here, on the other hand, we show that at least immediately 
outside their nuclear regions, Seyfert Is and 2s are very sim- 
ilar, and that dilution is really a nuclear effect present only 
in type Is. To further illustrate this, Fig. [9] repeats the com- 
parison performed in Fig. [T] but this time using off-nuclear 
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Figure 9. Similar to Fig.[T]but for ofT-nucIear measurements: (a) 
1 kpc data, (b) furtliest point for which we have signal (see text 
for details). Scale and areas as in FiglTI 



W^CaT measurements. Depending of the galaxy redshift, the 
furthest off-nuclear measurements correspond to less than 1 
kpc to ~ 3 kpc. In Fig.|9^ we proceed as follows: in the first 
cases (nearest galaxies) we take the furthest points from the 
nucleus (for instance, NGC 1433 in Fig.|8]), and in the second 
cases (furthest galaxies) , we firstly obtained the interpolated 
CaT values at exactly 1 kpc taking into account the data at 
each side of the nucleus, and then we calculate the average 
of the interpolated values (for instance, NGC 3081 in FiglH}. 
We see that the differences between Seyfert Is and 2s disap- 
pear when considering these off-nuclear values, which we call 
"<lkpd' data. Fig. |9j3 does the same thing, but using the 
values corresponding to the furthest points from the nucleus 
for which we have signal, and we call them "furthesf data. 
Numerical results of Wcs.t spatial measurements are sum- 
marized in Table [5] A KS test applied to the distributions 
yields probablitities of p = 0.46 ("<JA;pc") and p — 0.84 
{"furthesf). We see that even for short distances (< 1 kpc) 
there is no significant dilution of CaT in Seyfert Is, as we 
expected from the regions outside the nuclei. 

It is also interesting to point out that these results are 
essentially identical to those obtained from the spatial anal- 
ysis of the CaK strength by Cid Fernandes etal. (1998) and 
Storchi-Bergmann etal. (1998, their Fig. 1). Hence, the FC 
which dilutes the optical spectrum in Seyfert Is is the same 
which decreases the CaT strength in the near IR. This com- 
ponent is not present in type 2s. 

The behavior of WcaT{r) for Starburst + Seyfert 2 com- 
posites (Figs. |8^ and f), on the other hand, contrasts with 
that derived from Wca.K radial profiles. Contrary to the 
strong spatial dilution by young stars detected in the CaK 
and other optical lines, the CaT profiles show very little, 
if any, variation as a function of distance to the nucleus in 







Nuclear [A] 


<1 kpc [A] 


Furthest [A] 


Seyfert 1 


- min 


2.7 


4.3 


5.0 


Seyfert 1 


- mod 


4.8 


6.5 


6.6 


Seyfert 1 


- max 


6.7 


8.3 


9.0 


Seyfert 2 


- min 
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Table 2. Statistics for Spatial WcaT measurements: minimum, 
median and maximum values are shown for Seyfert 1 and Seyfert 
2 and for nuclear and off-nuclear data. For off-nuclear data, we 
considered < 1 kpc data, which refers to the values at 1 kpc or less 
from the nucleus, and furthest data, which are the corresponding 
values of furthest points from the nucleus of each galaxy for which 
we have signal. 



composite systems. In NGC 7582 there is a marginal hint of 
dilution. From optical work we know that this galaxy hosts 
a very dusty central Starburst, so it could represent a case 
where the Starburst light does contribute to the fiux in the 
CaT range with a non negligible fraction. However, even in 
this most favorable case the nuclear CaT is diluted by less 
than 15% with respect to its value at r ~ ±500 pc. In con- 
trast, the optical absorption lines are much more diluted: 
At the wavelengths of the CaK, G-band and MgIb lines, the 
dilution factors are 70, 60 and 57%, respectively. To bet- 
ter appreciate this strong difference, we invite the reader to 
compare Fig. [S^ with fig. 29 of Cid Fernandes etal. (1998), 
where the radial profiles of several optical absorption lines 
in NGC 7582 are shown. 

The difference in WcaT and Wca.K profiles is striking, 
but not surprising. In fact, our finding that the CaT does 
not mimic the CaK spatial behavior is fully consistent with 
the conclusion laid out in the section above that only under 
extreme conditions (very young age and large extinction) 
would a Starburst strongly dilute both CaK and CaT lines. 



5 KINEMATICS 

Our spectroscopic data (Paper I) provides estimates of two 
characteristic velocities: the typical velocity of stars through 
the host's bulge, cr*, and the typical NLR cloud velocity dis- 
persions, inferred from the width of the [SIII]A9069 line: We 
fitted the profile of [SIII]A9069 with one or two Gaussians, 
one representing the core of the emission line and the other, 
if necessary, representing the wings. The resulting FWHM 
of the core component divided by (81n2)^^'^ is what we call 
""[siii]- We have determined (T[siii] by this method for 31 
Seyfert and 4 Starburst nuclei. 

The distributions obtained for cr* and o"[siii] are shown 
in Fig. 1101 They show that both types of nuclei have essen- 
tially the same (statistical) values: for Seyfert Is cr* = 128 
± 37 km/s and crpni] = 111 ± 29 km/s, while for Seyfert 2s 
(J* = 134 ± 45 km/s and o-[siii] = 137 ± 76 km/s. Firstly, 
we will discuss the stellar velocity dispersions (cr*) and then 
we will analyze the NLR kinematics through (Jpnij and its 
possible link to cr*. 
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Figure 10. Distributions of stellar velocity dispersions (ct*, top) Figure 11. CaT values versus stellar velocity dispersions. Dashed 

and [SIII]-based gas velocity dispersions (crjgjjjj, bottom). Areas line represents Saglia's fit (see text). Symbols as in Fig. [2] 

as explained in Fig. [T] 



5.1 (T* versus Wca.T 

As reviewed in Section [T] recent work on normal galax- 
ies has dedicated much attention to the fact that WcaT is 
slightly anti-correlated with velocity dispersion all the way 
from dwarf to giant elliptical galaxies and bulges (Saglia 
etal. 2002; Cenarro etal. 2003; Falcon- Barroso etal. 2003; 
Michielsen etal. 2003). Even though our CaT data are not 
of the same quality as that employed in these studies, it is 
interesting to check whether Starburst and active galaxies 
follow the same WcaT-o"* as normal ones. 

Fig.[TT]shows the results of this test. Clearly, as a whole, 
our galaxies do not follow a well defined WcaT-o"* relation 
at all. The few normal galaxies in our sample (crosses), as 
well as the high-tj* Seyfert 2s, do tend to follow a similar 
anti-correlation found by Saglia et al (2002) for ellipticals. 
However, there is a large population of galaxies with low- 
er* and relatively low WchT which has no correspondance in 
previous studies of this relation. Excluding Mrk 3 and NGC 
7212 (which are diluted by EC, as discussed in l3.3p . the most 
deviant points are the Starbursts NGC 1140, NGC 3256 and 
NGC 7714, located in the bottom left part of the plot. Next, 
at WcaT between 5 and 6.5 A and cr* < 150 km/s, are 
the composite Starburst-Seyfert 2s nuclei, like NGC 5135, 
NGC 7582, IC 3639, Mark 1. It thus seems that the complex 
star formation histories of these objects are somehow related 
to their displacement with respect to the VKcaT-o"* relation 
defined for more well behaved systems. Interestingly, the 
normal galaxy NGC 205 is located away the anticorrelation 
stated for normal nuclei. This dwarf Elliptical galaxy has a 
low metallicity (Mateo, 1998), which indicates that it could 
also infiuence the CaT-tj* relation. 



5.2 Stellar vs. NLR Kinematics 

The question concerning the acceleration of gas in the NLR 
has been analyzed in several papers (NW, TDT) . As pointed 
out by Green & Ho (2005; GH), the NLR is smaU enough to 
be illuminated by the active nucleus and large enough to feel 
the gravitational forces of the bulge of the host galaxy. Com- 
parison between gas and stellar kinematics was previously 
carried out by NW using the [OIIIJA5007 line, who find that 
(T* and the EWHM of [OIII]A5007 are correlated, but with 
substantial scatter. Recent work on Seyfert Is (Jimenez- 
Benito et al. 2000) and SDSS galaxies (Heckman et al. 2004, 
GH) confirms this finding. The implication is that the NLR 
clouds are at least partly dominated by virial motions in the 
host galaxy's bulge. Given the existence of this relation, for- 
bidden line widths have sometimes been used as a surrogate 
for CT* in the absence of information on the stellar dynamics 
(Nelson 2000; Grupe & Mathur 2004). The case of Narrow- 
Line Seyfert 1 galaxies (NLSl) is particularly interesting. In 
a recent paper, Komossa and Xu (2007) have shown that, 
for these Seyfert galaxies, o-[oni] can be used as a proxy 
for CT* only after removing the (usually blue) wings of the 
[QUI] profile. In other words, cr[oiii] should correspond only 
to the core of the line, similarly to what we have done in the 
present paper in the fit of the [SHI] emission-line profiles to 
obtain o"[siii]. 

In Fig.[T2]we compare our values of Ufsiii] and cr* for all 
Seyfert nuclei. We note that two of them, namely Mrk 1 and 
NGC 1068, have o"[sni] much larger than cr*. Evidence of jets 
were found in these galaxies, which may explain the acellera- 
tion of the gas via shocks produced in the interaction within 
jets and NLR. Excluding these two nuclei, we have fitted 
all the other galaxies by linear relations given by Ordinary 
Least Square (OLS) regressions (Isobe et al, 1990). Solid 
line represents the OLS Bisector fit, while dashed lines are 
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Figure 12. Comparison of stellar velocity dispersions with our 
estimate of o-gas based on the width of the [SIII]A9069 line. Sym- 
bols as in Fig. [2] Ordinary Least Squares regressions are applied 
over all data but Mrk 1 and NGC 1068: solid, OLS Bisector 
fit; dashed, OLS{X\Y) and OLS{Y\X) (shallowest and steepest 
lines, respectively). 





All 


Sy 1 


Sy 2 


All (a + h) 


All (a) 


Slope 


0.93 


0.73 


1.03 


1.00 


1.00 


N 


29 


8 


21 


23 


20 


Rs 


0.37 


0.31 


0.40 


0.39 


0.30 


p 


0.046 


0.45 


0.07 


0.06 


0.20 



Table 3. Results of OLS Bisector fits for crjgjjjj-crt. Letters in 
columns 5 and 6 denotes quality flags (See text). 



OLS{X\Y) and OLS{Y\X) (shallowest and steepest lines, 
respectively). The OLS Bisector fit gives a slope of 0.93 (N = 
29, Rs = 0.37, P = 0.046). Taking into account just Seyfert 
2s, we find a slope of 1.03 (N = 21, Rs ~ 0.40, P = 0.07) 
and 0.73 for Seyfert Is, for which there is not a clear cor- 
relation [Rs ~ 0.31, P = 0.45). The slopes and correlation 
coefficients do not vary substantially when we exclude worst 
quality data. Using the a = very good, 6 — good, and c = 
not-too-bad quality flags adopted in Paper I, we find for 
a + & a slope of 1.00 (N = 23, Rs = 0.39, P = 0.06), and 
1.00 for a sources alone (N = 20, Rs = 0.3, P = 0.2). All 
these is summarized in Tabled Interestingly, there is also a 
tendency for 4 Starburst nuclei, with a slope of about 0.56 
and Rs ~ 0.81 (not shown here). 

We repeated this analysis using ""[om] = 
FWHM[oiii]/(81n2)^/^ taken from the literature (Whittle, 
1992; Bassani etal. 1999; Schmitt etal. 2003). NGC 1068, 
Mrk 1 and Mrk 78 present very different values of cTjoiii] 
and (T*. Further excluding the worst quality data ("flag c" 



AU Sy 1 Sy 2 

Slope 1.48 1.13 1.60 
N 28 12 16 

Rs 0.42 0.81 0.33 

Table 4. Results of OLS Bisector fits for cjoiiij-o"* 



in Whittle 1992 and our flag d) the OLS Bisector fit gives a 
slope of 1.48 {Rs = 0.42) for the whole sample (28 objects), 
while splitting into Seyfert Is and 2s the slopes are 1.13 
{Rs = 0.81) and 1.60 {Rs = 0.33) respectively. This is 
described in Table |31 We point out that for this analysis 
we used FWHM data of [OUT] line, which are (in median) 
larger than our estimates of o-gas through the FWHM of the 
core of [SHI] line, so that the slope of a fit applied to [OIII] 
data yields larger values than those obtained for [SIII] data. 
We conflrm this fltting o^^^-^^ vs (^^q^] , which gives a slope 
of ~ 0.85 {Rs ~ 0.44, 20 objects). 

There are, therefore, correlations between NLR and 
stellar motions, but with considerable scatter. One way to 
study this scatter is to look at the residuals Ojsjjjj /cr* about 
unity. Fig. [13] shows this distribution, which is a moder- 
ately peaked (kurtosis ~ -1) gaussian for the majority of 
Seyfert 2s. Without taking into account the two galaxies 
with abnormally high gaseous velocity dispersion (for which 
it was shown that there are several kinematical compo- 
nents), Cigni]/""* ~ 0-9''' =t 0-32 (average and rms disper- 
sion), while for Seyfert Is and 2s o-jgjjjj/tj* are 1.02 ± 0.40 
and 0.95 ± 0.29, respectively, i.e., /y^sm]/'^* same for 

both types of objects. Fig. [T3] is similar to fig. 7 of NW, 
which uses [OIII]-based gas kinematics. Although NW have 
twice as many objects, we note that our distribution is more 
concentrated around unity, mainly because we have FWHM 
determinations of gaussians fitted to the core of the [SIII] 
line. 

All this study shows that the [SIII] core emitting zone 
of the NLR is dynamically linked to the gravitational bulge 
potential, a result similar to that obtained by Komossa and 
Xu (2007) for NLSl galaxies. However, we see that there 
are possible secondary factors besides gravitational potential 
that could explain the scatter about unity in NLR kinemat- 
ics. Shocks and other acceleration mechanisms could account 
for the high velocity clouds observed in a few galaxies, but in 
general a complete study of gas-stellar kinematics remains 
to be done. 

A detailed analysis of stellar and gaseous kinematics of 
NLR was made by GH taking precisely cfjgjjjj /cr* as the cen- 
tral quantity of their study. They found that none but one of 
the nuclear activity indicators is correlated with i^fsiu]/'^*' 
only with log a*, which is indirectly a measure of bulge mass, 
there is a hint of relationship. Although we have few points, 
we note this trend in our data, which is shown in Fig. [14] 
together with GH's best flt. Again, applying OLS Bisector 
fit, we obtain log o-jgjjjj /ct* — — l.llogcr* -I- 2.2, with Rs — 
—0.55. A least squares fit gives log f j^^^jj /a* — —0.56 log a*-|- 
1.13 {Rs ~ -0.55, P = 0.022), nearly the same relation ob- 
tained by them. Only for Seyfert 2s there is a weaker (but 
real) trend, log cr^gj^j /ct* = -0.41 log cr* 4-0.81 (N = 21, 
Rs = —0.37, P — 0.09), while for Seyfert Is the correlation 
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Figure 13. Distribution of gas to stellar velocity dispersions ra- 
tio, (Tjgjjjj /(T^. Areas as explained in Fig. [T] 



0.4 




2 2.2 

log(CTjkm.s-i]) 

Figure 14. Deviation from o-isiu]/'^* with cr*. Solid line is the 
least square fit for all Seyferts. Dotted line, Green&Ho's fit. Sym- 
bols as in Fig. [2] 



is moderately strong, log CTj^jjjj /cr* = —0.74 log cr* + 1.50 (N 
= 8, i?s = 0.75, P = 0.036). 

We checked whether crfsm] correlates with other 
properties, like morphological and activity types, redshift, 
inclination, etc., but (like GH) found no strong trend. The 
absence of a trend with the inclination of the host galaxy 
suggests that the gas in the NLR has random velocities and 
so the velocity field is similar to that of the stars, as previ- 
ously found by NW (their Fig. 11). 

Our data thus support the existence of an anticorrela- 
tion between (T[siii] and a*, albeit with a large scatter. Using 
sigmttjgjjjj as a proxy of cr* tends to underestimate velocity 
dispersion in bulge of massive galaxies. In any case, given 
the existence of a correlation between both (gas and stellar) 
velocity dispersions, and taking into account that [SIII] is a 
high ionization line produced very near the active nucleus 
and often with strong wings like [OIII], our study supports 
the idea of using o-jg^j^ as a kinematic tracer of the NLR. 



5.3 Spatially resolved kinematics 

As an extension of our spatial analysis, we report results for 
the variations in kinematical properties of our sub-sample. 

Fig.[T5]shows examples of spatial variations in a*. Some 
recent works (Barbosa et al. (2006) , Garci'a-Lorenzo et al. 
1999; Emsellen etal. 2001; Marquez etal. 2003) have dis- 
cussed the possibility of a sigma-drop detection, i.e., that 
a* is appreciably lower towards nuclear regions due to the 
presence of a star-forming cold nuclear disc. We see a nu- 
clear drop in ct* in only two galaxies (Fig. 115b , and b). For 
some objects there is even a hint of a rise in cr^ in their cen- 
tral parts (e.g.. Fig. I15fc and d). We have to keep in mind, 
however, that our spectra do not have the spatial resolution 



needed to detect this behavior in cr* , which would require a 
resolution of a few tens parsecs. 

Another interesting feature that can be seen for some 
of our objects is a symmetric off'-nuclear drop in cr*, as in 
NGC 4593 (Fig.Iin}?). This is in agreement with the detailed 
integral field spectroscopy of this object made by Barbosa 
et al. (2006), who found a partial ring of low-cr* nearly at the 
same distance from the nucleus. A similar behavior is seen 
in NGC 6951 (Fig. I15F 1. U-band images of these galaxies 
(Mufioz Martin et al, 2007) shows the presence of young 
stars at the locations where the off nuclear a* drops are 
observed. 



6 SUMMARY 

In Paper I we presented a spectroscopic atlas of 78 galax- 
ies in the region of CaT lines, and measured the equivalent 
widths of CaT and [SIII]A9069 lines, and stellar and gaseous 
velocity dispersions. In this Paper we used these data, as 
other measurements, to investigate the CaT strength and 
the kinematical properties of Seyfert nuclei. 

Our conclusions can be summarized as follows: 

(i) The equivalent widths of CaT and [SIII] are diluted 
in Seyfert Is, due to the presence of a non-stellar compo- 
nent, while in Seyfert 2s there is no sign of dilution. We 
show that the nuclear dilution of CaT lines in Seyfert Is 
and the non observed dilution in Seyfert 2s holds spatially, 
i.e., within approximately the central kpc. At this distance 
the effects disappears, thus implying that the stellar lines 
are diluted due the nuclear non-stellar continuum directly 
seen in Seyfert Is. 

(ii) The CaT strength turns out not to be a good tracer 
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Figure 15. Examples of ct* spatial variations. Dotted lines mark 
—200 to +200 pc region; dashed lines, the —1 to +1 kpc region. 



for the mean ages and metallicities of metal rich stellar pop- 
ulations. Nonetheless, its combination with CaK line yields 
a useful constraint on the nature of the continuum emission 
from optical to near-IR wavelengths. 

(iii) We show that the location of Seyfert 2s in the CaT- 
CaK plane could be explained satisfactorily by stellar pop- 
ulation synthesis models, by considering mixtures of an old 
plus very young stellar populations with a extinction of just 
two or three magnitudes. The hypothesis of Seyfert 2 nuclei 
composition as a mixing of an old stellar population plus 
power-law central source cannot account for the majority of 
the data. 

(iv) There is a correlation between nuclear cr* and <T[siii] , 
as well as an anticorrelation between cr[siii]/a* and a*, both 
with substantial scatter. This means that care must be ex- 
ercised when using (Jpin] as a proxy for cr*. These results 
are compatible with previous results of NW and GH. 
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